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We calculate the thermal conductivity of electrons and muons K efl produced owing to electromag- 
netic interactions of charged particles in neutron star cores and show that these interactions are 
dominated by the exchange of transverse plasmons (via the Landau damping of these plasmons in 
nonsuperconducting matter and via a specific plasma screening in the presence of proton supercon- 
ductivity). For normal protons, the Landau damping strongly reduces re eM and makes it temperature 
independent. Proton superconductivity suppresses the reduction and restores the Fermi-liquid be- 
havior K eM oc T" 1 . Comparing with the thermal conductivity of neutrons k„, we obtain K eM > K n 
for T > 2 x 10 9 K in normal matter and for any T in superconducting matter with proton critical 
temperatures T cp > 3 x 10 9 K. The results are described by simple analytic formulae. 
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I. INTRODUCTION 

Neutron stars are compact stars of mass ~ 1.4 M and 
radius <~ 10—15 km, where M Q is the solar mass. A 
neutron star is thought to consist of a thin crust (~ 1% 
by mass) and a bulky core [E 0> @, 0- The crust is 
composed of a nonuniform matter containing atomic nu- 
clei. The core extends from the layer of the density 
p « 0.5 po to the stellar center [p ~ (10 — 20) po], where 
po ~ 2.8 x 10 14 g cm -3 is the mass density of saturated 
symmetric nuclear matter. It is widely thought that the 
core is built of neutron-rich, strongly degenerate, uni- 
form nuclear matter of supernuclear density. This makes 
neutron stars unique natural laboratories of such mat- 
ter. The composition and equation of state (EOS) of this 
matter are largely unknown but will hopefully be strongly 
constrained soon by comparing theoretical models of neu- 
tron star structure and evolution with observations. 

We will consider the thermal conductivity of neutron 
star cores. The thermal conductivity of neutron stars 
is needed for modeling the thermal structure and evo- 
lution of these stars and related phenomena [!, 0. In 
particular, it is important for cooling of isolated neutron 
stars and for thermal structure of transiently accreting 
neutron stars (in soft X-ray transients). Accreting stars 
can be heated by pycnonuclear reactions when accreted 
matter sinks into the deep neutron star crust under the 
weight of newly accreted material 0, H, • Thermal con- 
duction spreads this heat over the star - to the surface 
(resulting in surface photon emission) and to the core 
(producing neutrino emission from the core). In this way 
the observable surface thermal emission depends on the 
neutrino luminosity of the star, and, hence, on its in- 
ternal structure. The thermal conductivity of neutron 
star cores is especially important for studying cooling of 
young neutron stars (of age t < 100 years), where the in- 
ternal thermal relaxation is not yet achieved (e.g., Refs. 
[lol dH). The thermal conductivity can also regulate 
thermal relaxation of a neutron star in response to pul- 
sar glitches (see, e.g., [12j, and references therein). 



We will focus on models of neutron star cores composed 
of neutrons (n), with an admixture of protons (p), elec- 
trons (e), and muons (ji). Our results are easily general- 
ized for matter containing hyperons. In all these cases, 
the thermal conductivity n can be subdivided into two 
relatively independent parts [H, [3| , 



(1) 



where Kb is the conductivity of neutrons (and possibly 
other baryons) mediated by strong-interaction collisions 
of heat-conducting particles (see, e.g., [l5[, and references 
therein), while K efl is the conductivity of electrons and 
muons mediated by electromagnetic interactions. The 
conductivity K efx has been calculated in several papers, 
particularly, in Refs. [TEHEES- 



However, all previous calculations of K e/J in neutron 
star cores have neglected an important effect of the Lan- 
dau damping of electromagnetic interactions owing to the 
exchange of transverse plasmons. In the context of trans- 
port properties of dense matter this effect was studied 
by Heiselberg and Pethick for a degenerate quark 
plasma. Recently we have included the effect of the Lan- 
dau damping into the thermal conductivity of strongly 
degenerate electrons determined by electron-electron col- 
lisions with the electron plasma screening p^ |. Similar 
effects have also been studied by Jaikumar et al. [2(| for 
neutrino bremsstrahlung radiation via electron-electron 
collisions in neutron star crusts and cores. Here we re- 
consider K efl in a neutron star core including the effects of 
the Landau damping and of the specific plasma screening 
induced by proton superconductivity. 



II. ELECTRON-MUON THERMAL 
CONDUCTIVITY IN NORMAL MATTER 

The electrons and muons in a neutron star core con- 
stitute strongly degenerate almost ideal Fermi gases 0. 
The electrons are ultrarelativistic. The muons can be 
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absent in the outermost part of the core; they are non- 
relativistic just after their creation threshold and become 
relativistic at higher densities. The neutrons and protons 
constitute strongly interacting Fermi liquids. The neu- 
trons are mainly nonrelativistic (although they become 
mildly relativistic near the center of a massive neutron 
star); the less abundant protons are typically nonrela- 
tivistic. The neutrons and protons can be in superfluid 
states; superfluidity of protons means their superconduc- 
tivity. In this section we consider nonsuperfluid (normal) 
matter. We analyze the effects of proton superconduc- 
tivity in Sec. IIIII We closely follow the derivation of K efJl 
by Gnedin and Yakovlev [171 ] (but introduce the Landau 
damping); thus we omit technical details. 

The thermal conductivity K efJi is calculated from a 
coupled system of linearized Boltzmann equations for 
the electron and muon distribution functions f e and / M , 
which we denote collectively by f c , with c = e or /j. 
These distributions slightly deviate from the equilibrium 
Fermi-Dirac distributions f c owing to the presence of a 
weak temperature gradient VT, 



fc = /c-* 

f c = \ exp 



dfc 
d£ c ' 

£ c - Vc 
k B T 



1 



(2) 



where e c is the particle energy, p, c is its chemical poten- 
tial (the rest-mass term included), T is the temperature, 
ks the Boltzmann constant, and $ c measures a devia- 
tion from equilibrium. For calculating n efJi , the electrons 
and muons are treated as the only heat carriers which 
experience electromagnetic interactions between them- 
selves and with any charged baryons in dense matter 
(with protons, in our case). The charged baryons are 
assumed to obey equilibrium Fermi-Dirac distributions 
(see, e.g., Ref. [IT] ) . Nonequilibrium parts of the elec- 
tron and muon distributions are found using the standard 
variational approach with the simplest trial function, 



*c = -T C (E C - fie) v c ■ VT/T, 



(3) 



where v c is the velocity of particles c and r c is their ef- 
fective relaxation time. 

The resulting electron and muon thermal conductivity 
can be written in a familiar form 



= K 
IT 



2 k 2 Tn P i 



3m* 



3m* 



(4) 



where n e and are the partial thermal conductivities of 
electrons and muons, respectively; n e and are number 
densities of these particles; m* = fi e /c 2 and m* = /Lt^/c 
are their effective masses. In neutron star cores at beta 
equilibrium one has \i e = /i M and m* = m* pj. Note 
that the thermal conductivities n e and presented by 
Gnedin and Yakovlev [13 [their Eqs. (61) and (65)] con- 
tain also an additional factor C w 1.2 which brings the 



employed simplest variational solution closer to the exact 
solution. We do not introduce a similar correction here. 
Its calculation is more complicated than in Ref. [13], be- 
cause now we include the Landau damping of transverse 
plasmons. The Landau damping introduces non-Fermi- 
liquid behavior into the thermal conductivity (Sec. Ill Bj) . 
while the standard technique for calculating the correc- 
tions is developed for traditional Fermi liquids. In Sec. 
Ill CI we argue that the Landau damping of transverse 
plasmons reduces the difference between the exact and 
simplest variational solutions. Thus, neglecting the cor- 
rection factor C in the expression for can introduce 
uncertainties <20% which are quite acceptable for the 
thermal conduction problem in neutron star cores. 

The effective electron and muon relaxation times can 
be written as 1171 



VcVv 



V e Vn 



(5) 



where 



v e = yv, 

i 



fj,e i u 'fj,p 



(6) 



are the total effective collision frequencies of electrons 
and muons with all charged particles i (which are elec- 
trons, muons and protons, in our case); v e i and v^i are 
familiar partial effective collision frequencies (derived in 
a standard manner from respective linearized collision 
integrals), while v' efl and v'^ e are two additional effec- 
tive collision frequencies 17] which couple heat transport 
of the electrons and muons. All these collision frequen- 
cies can be expressed as multi-dimensional integrals over 
momenta of colliding particles. Since all particles are 
strongly degenerate, their momenta can be placed at ap- 
propriate Fermi surfaces whenever possible. In this way, 
for collisions of nonidentical particles 1 and 2 we obtain 
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while for collisions between identical heat carriers we 



have v cc = v\i 



Here, vp\ and ppi are, respec- 



tively, the Fermi velocity and the Fermi momentum of 
particles 1; 



xM/(l,2|l',2')/i/ 2 (l-/0(l-/2) 

x (v a ■ v f3 ) (e a - fj, a ) (ep - up) (8) 

is the normalized collision frequency; 

H/(l,2|l',2') = 4(2^) 4 5(e 1 +e 2 - £ ' 1 -4) 

y.5{ Pl +p 2 -p l 1 -p' 2 )\M l2 \ 2 (9) 
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is the differential transition probability for a scattering 
process P1P2 v'\P'%\ I-W12 1 2 is the squared matrix ele- 
ment summed over particle spin states (it includes also 
the symmetry factor to avoid double counting of the same 
collisions events) ; primes refer to particles after the colli- 
sion. In Eqs. (J?])-© we use the system of units in which 
c = = h = 1. The same system will be used below 
unless the contrary is indicated. 
In the absence of muons we have 



(10) 



(q <C PFi, u <C VFiPFi) one has (e.g., Alexandrov, Bog- 
dankevich and Rukhadze (22J) 



i 

n* = £ 



TfVFi 

— Xt{xi), 



(12) 



where Xi = u>/(qvFi), a = e 2 /hc ~ 1/137 is the fine 
structure constant, and 



Dynamical screening of electromagnetic 
interaction 



The physics of the dynamical plasma scree ning is thor- 
oughly analyzed by Heiselberg and Pethick [lq . These 
authors consider quark-quark collisions in a quark plasma 
through one-gluon exchange in the weak-coupling limit. 
Such collisions are similar to electromagnetic scattering 
of charged particles in an ordinary plasma. Electromag- 
netic interactions of muons and electrons in neutron star 
cores are usually accompanied by small momentum and 
energy transfers which greatly simplifies the theory. The 
squared matrix element for small energy transfers in a 
collision P1P2 — * v'\P'i °f nonidentical particles 1 and 2 



, , x 2 x(l-x 2 ) . fx + l\ . „. 

»(-) = y + ^-^ ln (;^-Tj- ( 13 ) 

For typical conditions of very strong degeneracy in neu- 
tron star cores, it is sufficient to use the expressions for 
II; and Il t in the limit of q <C PFi and x% <C 1, in which 
Xi(x) ~ 1 and Xt( x ) ~ iirx/A. In this limit, 
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where hq = p[ p\ is a momentum transfer, fku = 
e'i — £1 is an energy transfer (in standard physical units), 



J, 



(4 



(0) 



Ji'i) — (u' i J 1 ' u i) is the transition 4-current 
(u=0,. . . ,3), 7" is a Dirac matrix, is a normalized 
bispinor (ttjHj = 2to^), Ui is a Dirac conjugate (see, e.g., 
Berestetskii, Lifshitz and Pitaevskii [H]]); Jun is the com- 
ponent of Ji'i transverse to q. The longitudinal com- 
ponent of Ji'i (parallel to q) is related to the time-like 

(charge density) component J$ via charge conservation. 
It is excluded from Eq. (fTTj) with the aid of the continuity 
equation as explained by Heiselberg and Pethick [18[ . For 
collisions of identical particles (ee and fj,fj, in our case) , the 

matrix element contains two parts, M\% — + , 
which correspond to two channels, (1 — > 1'; 2 — > 2') and 
(1 — > 2'; 1' — > 2). However, in the small-momentum- 
transfer approximation, the interference term is small, 
both channels give equal contributions, and the relation- 
ship (fTTj) is not violated. 

Equation (fTTj) contains the polarization functions III 
and n t which depend on u> and q and describe plasma 
screening of interparticle interaction by longitudinal 
and transverse plasma perturbations (plasmons), respec- 
tively. In the random phase approximation (RPA), these 
functions are the sums of terms for all charged particles 
i (electrons, muons and protons). In the classical limit 



where uJi — (47re 2 rii/m*) 1 / 2 is the plasma frequency 
of particles i, and m* is the effective particle mass at 
the appropriate Fermi surface. We have already defined 
m* and m* for electrons and muons which form almost 
ideal Fermi gases. Their effective masses differ from the 
bare masses owing to relativistic effects. As for protons 
(i = p), which are nonrelativistic in neutron star cores, 
their effective mass m* differs from the bare proton mass 
due to strong interactions with surrounding nucleons. In 
the approximation (fl"4"j) and (fT5|) one can also neglect u> 2 
in the denominator of the second term in (JTTJ) . Equations 
(fT"3"j) are strictly valid for Fermi gases, which is a good ap- 
proximation for the electrons and muons, but the protons 
constitute a Fermi liquid. Fortunately, the asymptotic 
expressions (fl4|) and (fT5|) . sufficient for our calculations, 
remain valid for the Fermi liquid. This is because a ki- 
netic equation for quasiparticlcs in the Fermi liquid at 
low ui is similar to a familiar kinetic equation for Fermi 
gases (see, e.g., Ref. (23|). 

In Eqs. (fi"4"|) and (fT"5"|) we have introduced qi and qt 
[cm -1 ] defined as (in standard physical units) 



QI 



4a ^ crnlppi 

7T / 



£ 

i 

£■ 



h 2 



h 2 



4a v^lfo 

h 2 



£ 



(16) 
(17) 



Generally, we have qt < qi- If all charged particles were 
ultrarelativistic then q t — > qi- 
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According to Eqs. (fl"4"|) and (|15p . the screening func- 
tions for longitudinal and transverse plasmons are differ- 
ent. The screening via longitudinal plasmons (fl4|) is de- 
scribed by a real (nondissipative) polarization function 
Ilf = qf which is equivalent to a Debye-like (Thomas- 
Fermi) plasma screening with the screening length 1/qi- 
The screening of transverse currents is described by a 
purely imaginary (dissipative) function (|15p , correspond- 
ing to the Landau damping via collisionless absorption of 
transverse plasmons by plasma particles. Calculating the 
electron and muon thermal conductivity in neutron star 
cores, the authors of Refs. [HI, [13] have neglected the 
Landau damping and erroneously used the same static 
Debye-like screening, IT t = qf, for transverse currents as 
for longitudinal ones. We will show that this approxima- 
tion strongly overestimates the thermal conductivity of 
electrons and muons. 

The squared matrix element (fTTjl for collisions between 
heat carriers c and scatterers i can be written as 



tides to 



\MmV = 



16n 2 a 2 
m*. m$ 



<P = W\\ + ¥>-l + ¥>_l|| 



(18) 



where tp is a dimensionless function determined by the 
exchange of longitudinal plasmons transverse plas- 

mons (<p_i_), and by both interactions (ip\\±). The leading 
terms in the weak screening approximation are 



n = 
<p± = 



*2 *2 

m* ml 
(q 2 +q?) 2 ' 

PFcPFiQ 2 COs2 t 
q e + A 6 

2m*m*pFiPFcq 4 cos < 

(q 2 + qf)(q 6 + A 6 ) 



(19) 
(20) 
(21) 



where A = (iTLoq 2 /4) 1 / 3 and 4> is the angle between pi and 
P2- For collisions of identical particles (i — c), the func- 
tion <p is generally different but the difference is negligible 
in the leading-order weak screening approximation. 



B. Calculation of collision frequencies 

Calculating the integrals (JSJ) for a strongly degenerate 
matter, we can use the standard angular-energy decom- 
position (e.g., Ref. [1]). It consists in placing momenta 
of colliding particles at respective Fermi surfaces, when- 
ever possible, and in decoupling energy and angular in- 
tegrations using d 3 p = m*pp de dfl. The squared ma- 
trix element in ^ depends only on relative orientations 
of momenta of colliding particles which greatly simplifies 
angular integrations. Performing similar transformations 
as in Gnedin and Yakovlev [17| we can reduce Eqs. ([7]) 
and ([5]) for the collision frequencies of nonidentical par- 



12a 2 



Tt^T 3 pFcTn 

!2a 2 p F 



tt^pI^ 



t \ JciX Jcil'ji 
(JcQ ~ Jci2')- 



(22) 
(23) 



Here. 



Jou = I I I I d£i d £ ; de 2 de' 2 hf 2 (l - f[){\ - f 2 ) 
x(ei - m)(ek - Mfc)^nfc(w) (24) 

is an energy integral and Ink(^>) is a remaining angular 
integral, which depends on a collision energy transfer u 
owing to the dynamical screening in Eq. (fT5|) . The effec- 
tive collision frequency v cc between identical heat carriers 
can be calculated as v cc = v C i + v' ci with i = c. 
In the weak screening approximation we obtain 



Iqi = Inv = dq / defxp, (25) 
Jo Jo 

In2 = Iq.2' = I dq d(j) (p cos 0, (26) 
Jo Jo 

with q m = mm.{2pFi,'2pFc}- The integration over (f> is 
then trivial. From Eqs. (fl"9|) - (j2Tj) we see that Im consists 
of two terms associated with the longitudinal and trans- 
verse plasmon exchange, respectively, while Iq,2 contains 
only the mixed term, 



1 m 



Ifl2 — I i 



U2 



(27) 



This means, that v C i is determined by the contributions 
from longitudinal and transverse screening, but v' ci is de- 
termined by the mixed term. 

We have calculated the angular integrals in the weak 
screening approximation, retaining the leading terms in 
series expansions over qi or A. The result is 



l Q2 



KPFiPFc 1 

3 k 



2 ' 

qf 



2 *2 *2 

7T ml TfU 



4qf 



2 

7T 



1/3 



PFiPFc m l m c 

qf 9t 2/3 M 1/3 



(28) 
(29) 
(30) 



Using the standard technique (e.g., [T3|) we can reduce 
the energy integration to 



Jcil j cil' 



Jc.i.2 — Jcil' — — 



Iqi(ui)uj 3 exp(— lo/T) duj 
[1 - exp(-^/T)] 2 ' 
I n2 (u)uj 3 exp(-uj/T)duj 
[1 - exp(-^/T)] 2 



(31) 



(32) 
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The final oj integration gives analytical expressions for 
the collision frequencies (in standard physical units) 



,11 



24C(3) a 2 k B T p 2 FtPFc 



7r 3 hem* h 2 q 2 ' 
4ir 2 a 2 k 2 B T 2 c A mfmf 

5 hc 2 m* p Fc h 3 qf 
12£ a 2 {k B Tf' 3 cm* c p 2 Fl 



7T 3 h(hcq t ) 2 / 3 h 2 p Fc q 2 



18.52, 



(33) 
(34) 

(35) 

(36) 

(37) 



where £(z) is the Riemann zeta function, and T(z) is the 
gamma function. Equations ([I])-© and (f33"|) - (f36"|) fully 
determine the electron and muon thermal conductivity 
of a nonsuperconducting matter calculated in the weak 
screening approximation. 

Owing to a strong degeneracy, for typical conditions in 
neutron star cores 



< uL < vl 



(38) 



Because v' and v' are negligible, electron and muon 
heat transports (c = e and /i) become formally decou- 
pled. In this approximation from Eq. ([5|) we obtain 



6C(3) ak B Tp Fc 
7r 2 hem* 



(39) 



Here the electron and muon thermal relaxation times r e 
and are mostly determined by the exchange of trans- 
verse plasmons (by the Landau damping) and are ex- 
pressed only through the parameters of heat carriers (e 
or fi). 

The thermal conductivity of the electrons or muons 
becomes temperature independent, 



7T 2 k B C P 2 Fc 



54((3) h 2 c 



(40) 



being solely determined by the Fermi momentum p Fc of 
heat carriers. 

The remarkable simplicity of the thermal conductiv- 
ity of strongly degenerate particles in this regime, where 
electromagnetic interactions are dominated by the Lan- 
dau damping of transverse plasmons, was pointed out 
by Pethick and Heiselberg [l8| for a plasma of mass- 
less (ultrarelativistic) quarks. The authors stressed that 
the temperature dependence of this thermal conductiv- 
ity strongly deviates from the dependence k oc T^ 1 in 
an ordinary Fermi liquid, because the plasma screening 
due to the exchange of transverse plasmons depends on 
energy transfer hio in collision events. Were particle col- 
lisions dominated by the exchange of longitudinal plas- 
mons [y — i/"), the thermal conductivity would show 
the traditional Fermi-liquid behavior (as was obtained in 



Refs. [H, [H, [13] , whose authors erroneously assumed the 
Debye-like plasma screening, the same as determined by 
the exchange of longitudinal plasmons). 

Recently we [l9| have calculated v eei taking proper ac- 
count of the electron plasma screening via the exchange 
of longitudinal and transverse plasmons, for degenerate 
electrons of any degree of relativity. We have used more 
general electron polarization functions given by Eqs. (fT2|) 
and (fT3")l [rather than by the asymptotic limits (fT4")l and 
(|15p employed in the present paper]. We have applied 
our results for studying the electron thermal conductiv- 
ity in a neutron star crust, but after adding the plasma 
screening by muons and protons [see Eqs. (fTtJ)) and (fTT)) ] 
they become valid for nonsuperconducting neutron star 
cores. Replacing in v ee electrons by muons we can obtain 
the expression for v^. 

Our present results include electromagnetic interac- 
tions of electrons, muons and protons (not only v ee or v, t u , 
which could be extracted from the results of Ref. [l9(); 
they are valid for degenerate particles of any degree of rel- 
ativity. Because the exchange of transverse plasmons is 
efficient for relativistic particles, it becomes unimportant 
if all charged particles are nonrelativistic. The presence 
of ultrarelativistic electrons in neutron star cores ensures 
the importance of the Landau damping there. 



C. Comparison with exact solution 

So far we have used a simplest variational thermal con- 
ductivity based on the approximation ([3]) with r c inde- 
pendent of e c . However, the actual energy dependence of 
the function $ c in the exact solution is more complicated. 
As already mentioned above, Gnedin and Yakovlev [ItJ 
introduced a factor C which corrects the variational 
thermal conductivity, K cxac t = CK var . For the electro- 
magnetic interaction wholly (and erroneously) attributed 
to the exchange of longitudinal plasmons, in the weak- 
screening approximation they obtained C ~ 1.2—1.3 and 
suggested to use C = 1.2. Their calculation was based on 
the exact transport theory in Fermi liquids developed by 
Sykes and Brooker [2~ij for one-component systems and 
extended for multi-component systems by Flowers and 
Itoh [l3[ and Anderson et al. (25|. Notice that if Gnedin 
and Yakovlev correctly calculated the thermal conduc- 
tivity solely determined by the exchange of longitudinal 
plasmons they would obtain the parameter A from their 
Eq. (62) equal to 1 and C = 1.2. The exact theory 0] 
cannot be directly applied to our case because the Lan- 
dau damping of transverse plasmons makes the matrix 
element ^-dependent, that is not characteristic for ordi- 
nary Fermi liquids. 

In this section we estimate C taking into account only 
the exchange of transverse plasmons in the weak screen- 
ing approximation for one type of heat carriers (recall 
that the electron and muon heat transports become de- 
coupled). We use the st raig htforward generalization of 
the exact theory [H, [H, HI |2f| . Instead of Eq. © we 
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write 



10 



4>, 



-Teff^(x) v c • VT, 



(41) 



where r e g is an effective relaxation time, while 'I'(x) is 
an unknown function of x = (e — /i c )/T. Following Ref. 
[26| we substitute (|41~j) into the linearized kinetic equation 
and obtain the equation for ^(x), 



1 + exp(x) 



dx sign(xi - x) ^(xi) - g(x) 
l + exp(xi) 1 — exp(x — xi)' 



with r e ff = 3/ (avFcT), and the thermal conductivity cal- 
culated as 



n c T c g 
ml 



dxx*(x)/ c (l-/ c ). (43) 



We have solved Eq. (|42| numerically, calculated K cxa ct, 
and compared it with the variational result (that corre- 
sponds to *var(x) = 7r 2 x/[18C(3)]). The difference be- 
tween the exact and variational results appears to be 
negligible, C = Kexact/^var ~ 1- In a more general case 
the electron (muon) scattering is determined by the ex- 
change of both transverse and longitudinal plasmons and 
we can expect that 1 < C < 1.2. 
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Figure 1: (Color online) Real and imaginary parts (thin and 
thick lines, respectively) of the proton transverse polariza- 
tion function Il t versus hui at small q in superconducting 
matter for the proton critical temperature T cp = 3 x 10 9 K 
and four values of the temperature, T = 10 s K (solid lines), 
2 x 10 9 K (dashed lines), 2.9 x 10 9 K (dotted lines), and 10 10 
K (dot-dashed lines). The plotted functions are divided by 
the imaginary part of life' in normal matter. 



III. EFFECTS OF PROTON 
SUPERCONDUCTIVITY 



Dielectric function 



In the previous sections we have considered nonsuper- 
conducting protons. However, it is well known that the 
protons in neutron star cores can be in a superfluid (su- 
perconducting) state (e.g., Lombardo and Schulze [271]). 
Let us include the effects of proton superconductivity. 

Microscopically, the onset of superconductivity after 
the temperature T falls below a (density dependent) crit- 
ical temperature T cp manifests itself in the appearance of 
a gap A = A(T) in the proton energy spectrum in the 
vicinity of the proton Fermi level. Instead of protons, it 
is now necessary to introduce Bogoliubov quasiparticles 
with energies (with respect to the Fermi level fi p ), 

e - n P = sign(p - p Fp ) A 2 + v 2 Fp (p - p Fp ) 2 . (44) 

It is generally believed, that proton pairing in neutron 
star cores occurs in the 1 5'o channel. In this case, the 
gap is isotropic (independent of orientation of proton 
momentum with respect to a quantization axis). For the 
Bardeen-Cooper-Schrieffer (BCS) pairing model the tem- 
perature dependence of A at < T < T cp can be fitted 
as [HI 

,.£.^(^-0^), (45) 

where t — T/T cp . 



Proton superconductivity affects electrodynamical 
properties of dense matter. The longitudinal polariza- 
tion function 11/ undergoes no significant changes (see, 
e.g., Ref. [29|]; a strong effect of superconductivity on 
longitudinal plasma screening described by the function 
Eq. (47) of Ref. [H is in error; one should set 
Z = 1 in equations of Ref. [I3|)- I n contrast, the 
(p) 

proton component HI of the total transverse function 
n t = n| e) + n| M) + l4 p) is considerably modified. It 
should be stressed that we discuss the electromagnetic 
polarization functions to be distinguished from polariza- 
tion functions associated with strong interactions (e.g.. 
Ref. [H). 

For our purpose, we employ the model of a Fermi gas 
of protons with BCS superconductivity in the so called 
Pippard limit appropriate to neutron star cores (huj <ti 
PFpVFp, hq < PF P , but £q > 1, with £ ~ hvF P /k B T cp 
being the coherence length). The polarization functions 
appropriate to this case are well known. They were cal- 
culated in the classical paper by Mattis and Bardeen [3l[ . 
In the limit of small u> and q, we are interested in, these 
results are gauge invariant [29(. The effects of strong 
(nuclear) interactions in the proton Fermi liquid (particu- 
larly, dragging of neutrons by protons) on the electromag- 
netic polarization functions are complicated and almost 
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not explored; they are beyond the scope of the present 
paper. According to M. E. Gusakov (private communica- 
tion) these effects should not greatly modify our results. 
Similar effects of strong interactions on the polarization 
functions of supcrfluid neutrons [301 ] lead to renormaliza- 
tion of vertices and take into account specific collective 
supernuid excitations. These effects can be pronounced 
for those q and w which are typical of such excitations 
(uj ~ <7«f), and the effects are small outside such q — uj 
domains. We expect that in our case (lu <C qvp) these 
effects are small and can be neglected. 

Therefore, using the classical polarization functions 
[3ll ] seems to be a good starting approximation. It gives 



4 q 



(46) 



where Q is a complex function of the normalized fre- 
quency uj — hid /A and the normalized gap y (or, alter- 
natively, of hw/ksT and y). The real and imaginary 
parts of Q are given by 



3?Q = 



= 



Q+l 

7T J tanh(|^ 

max [uj — 1,1] 

[E(E - Co) + 1] AE 

OO 

7T J jtanh (| (E + u)fj - tanh (| E^j } 



(47) 



[E(E- 



l]dE 



^{E 2 -l){{E + Co) 2 -l) 

cD-1 

7r6(a;-2) J tanh^l^ 
1 

x [E(u,-E)-l]dE 

y/ (E 2 — 1) ((d) — E) 2 — 1) ' 



(48) 



where 0(x) = at x < and O(x) = 1 otherwise. Here 
and below we again return to units c = h = k B = 1 
unless the contrary is indicated. 

In the nonsuperconding limit, y — > 0, only the imagi- 
nary part of Q survives, giving Q — iiruj and reproducing 
the asymptote (|15p determined by the Landau damping 
(dissipative plasma screening). In the opposite limit of 
strong superconductivity, y — > 00, the imaginary part of 
Q is small (because the gap in the proton energy spec- 
trum becomes much larger than characteristic values of 
hui ~ k B T). Then only nondissipative plasma screening 
is available (with Q = ir 2 ). These results are illustrated 

in Fig. [1] which shows the frequency dependence of the 

(p) 

real and imaginary parts of ILj in superconducting mat- 
ter. The presented functions are normalized by the imag- 
inary (Landau damping) part of Tl^ in normal matter. 
The critical temperature of proton superconductivity is 
taken to be T cp = 3 x 10 9 K {k B T cp = 259 keV). The 



frequency dependence is presented for four values of the 
temperature T. The highest temperature T = 10 10 K 
(k B T = 862 keV) refers to nonsuperconducting matter 
(y = 0, A = 0). The proton polarization is provided by 
the Landau damping. The displayed normalized polar- 
ization function has the imaginary Landau-damping part 
(the thick dot-dashed line which is the reference line of 
the figure); the real part is essentially zero. The dot- 
ted, dash-dot, and solid lines refer to progressively lower 
T = 2.9 x 10 9 , 2 x 10 9 , and 10 s K (k B T = 250, 172, and 
8.62 keV) in superconducting matter (with y = 0.570, 
2.26, and 52.6; and A = 142, 389, and 453 keV, respec- 
tively). One can observe the growth of the real part 
of the normalized polarization function, which becomes 
larger than the imaginary part, leading to a nondissipa- 
tive transverse plasma screening for T <C T cp . 

The total transverse polarization function of the 
plasma containing electrons, muons and superconduct- 
ing protons has the form 



IT = 



7TUJ ( r, A 

— q% — 5RQ 
4q I 7tuj 



+ i 



Qte 



It 



2 L <& — $K2 

^ TTUJ 



(49) 



where q 2 t — Aap 2 Fi /-K. In the case of strong supercon- 
ductivity, the main contribution to transverse screening 
comes from superconducting protons. When T decreases 
from T cp to zero, the transverse plasma screening turns 
from dissipative to nondissipative one. 

Thus, proton superconductivity modifies the trans- 
verse polarization function and the screening functions 
([20]) and (|2ip . One can show that in the weak screening 
approximation the angular integrals in superconducting 
matter can be presented as 



Im = lkiF\u/T,y,r), 
= In!F x Hu/T,y,r), 



(50) 
(51) 



where 1^ and J^J are the angular integrals in the non- 
superconducting matter, while 



F" 



ir\u\{r + 1) 

■nujr + A3Q| 

2 A3£Q 

1 arctan 

7r \Trujr + A3Q| 

[7r|w|(r + 1)] 1/3 \(nwr + A3Q) 2 



(52) 

1/3 



V3' 



arctan ■ 



\irujr -\ 
\nLur 



A3Q| 
f A9Q| 



AlkQ 
[7r|w|(r + 1)] 1/3 AMQ 



\TTU1T - 

2 



A3Q| (iruir + A3Q) + (AWtQ)* 



V3' 



1 Vkujt + A3QI 

- arctan — 

3 A5RQ 



1/6 



(53) 
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are the factors which take into account the effects of su- 
perconductivity. The parameter r is defined as 
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P 2 Fe 



Pf, 



Qtp 



PFp 



(54) 



It is a slowly variable function of the density in a neutron 
star core, constrained by the condition of plasma neutral- 



ity, n e 



n p . We have r 



1 in the absence of muons, 



and r > 1 if the muons are present. The maximum value 
r max = 2 1 / 3 w 1.26 is reached in the high-density limit 
in which the muons become ultrarelativistic. 

In the limiting case of strong superconductivity (y » 

1) 



F = 



2M(r + l) 

7T 2 A 



V3 V ^ 



1/3 



(55) 
(56) 



The collision frequencies involving nonsuperconduct- 
ing particles (ee, /i/i, and e/z) can be calculated directly 
from Eqs. (f3"Tj) and (j3"2"j) . taking into account the factors 
F 1 - and F^. It is convenient to present these collision 
frequencies in the form 



v cl R e (y,r), 



(57) 
(58) 



where and v' ci correspond to nonsuperconducting 
matter, while 



Ri(y,r) 



R x (y,r) 



l 



6C(3) Jo 
1 



■ 3 exp(-w/T) F ± duj 
[1 - exp(-u;/T)] 2 ' 



r(14/3)C(H/3) 



11 /3 exp (_ ti ,/T)Fll- L dw 



(59) 

[l-exp(- W /T)r 

describe the effect of proton superconductivity. In the 
case of strong superconductivity (i/ > 1) we obtain the 
asymptotes 

4?r 2 (r + 1) 



Rf{y,r) 



i?f (y,r) 



45C(3)y 



8tt 4 



r(14/3)C(H/3)15V3 



Try 



1/3 



.(60) 



We have calculated R'l (y, r) on a dense grid of y and 
r and fitted the results by the function 



( r + i)i/3 



[(r + l) 2 - 0.757y + (0.50651y) 



2ll/6 * 



(61) 



The fit reproduces numerical data, the above asymptote 
at y — > oo, and the evident condition i2|" L (0, r) = 1. The 
maximum relative fit error ss 8% occurs at y = 2.2667 
and r = 1.3. We do not present any separate fit formula 
for R^~(y, r) but obtain a practical expression for v^r s in 
the next section. 



B. Superconducting scatterers 

Aside of changing plasma screening, proton supercon- 
ductivity directly affects the partial collision frequencies 
v ep and v w . Now these collision frequencies describe in- 
teraction of electrons or muons with proton Bogoliubov 
quasiparticles. Because the number of such quasiparti- 
cles is not necessarily conserved, the collision frequencies 
should include contributions of three processes, (1) scat- 
tering (1,2 -► l',2'); (2) decay (1 -> l',-2,2'); and (3) 
coalescence (1,2,-2' — > 1'). Here, 1 and 1' refer to an 
electron or muon; 2, 2', —2, and —2' refer to Bogoliubov 
proton quasiparticles; —2 (or —2') stands for a proton 
quasiparticle with momentum and spin directions oppo- 
site to those for 2 (or 2'). A linearized collision integral 
I cp becomes 



7-coal 



L cp 



(62) 



where 



rsc 

cp ~ 



rdec 

cp 



rectal 

cp 



J dr/ 1 / 2 (i-/ 1 0(i-/ 2 

xl¥ c s p c (l,2|l',2')($i'-$i), 

\ J dT/ 1 (l-/_ 2 )(l-/ 1 ,)(l-/2 

x^ c d p oc (l|l',-2,2') 

^ / dT/i/a/.a^l-AO 



xH/ c c ° al (l,2,-2'|l') ($i 



(63) 
(64) 
(65) 



correspond to scattering (sc), decay (dec), and coales- 
cence (coal), respectively; dr = (27rh)~ 9 d 3 p2 d 3 pi> d 3 p 2 '- 
The factor | excludes double counting of the same colli- 
sion events. 

Let us focus on the differential transition prob- 
abilities W^(l,2|l',2'), W c d p oc (l|l',-2,2'), and 
W^p (1,2, — 2'|1'). In the case of scattering, according 
to Fermi Golden Rule, 



Wll = 2-K§{£x + e 2 -ev-e v 
where 



spins 



(l',2'|F|l,2) , (66) 



V = a J j dr 1 dr 2 [p c (ri)p p (r 2 )D\\(ri - r 2 ) 

-(J c (n)- J p (r 2 )) Dx(ri-ra)] (67) 

is the interaction operator, D\\ and D± being correspond- 
ing propagators. In the first approximation, protons in 
neutron star cores can be treated as nonrelativistic parti- 
cles. Then the second-quantized operators of the proton 
density and proton current are 



2m? 



(68) 
(69) 



9 



The psi-operator of proton field is 

$ = ^2 Xa exp(ip ■ r) a pa , 



(70) 



where a pa is a proton (not a Bogoliubov quasi-particle) 
annihilation operator and Xa 1S a un it basic spinor. Bo- 
goliubov quasiparticles are introduced through quasipar- 
ticle annihilation operators b pa basing on the Bogoliubov 
transformation 



b pa = u p a pa - sign(cr) v p a f _ p _ a , 
a pa = Upbpa + sign(cr)v p 6L f ,_ ., 



(71) 



where u p and v p are coherence factors which can be cho- 
sen in different ways. We will use the following set of 
these factors 



1 

sign(x) 
V2 



x 

1 + -, 

z 



z 



(72) 



where x = vf p {p — PFp)/T and z = (e — n P )/T. In 
the nonsuperconducting limit we have u p (y = 0) = 1 and 
v p (y = 0) = 0. With this choice the dimensionless energy 
spectrum of quasiparticles is 



T 



= sign(a;) \Jx 2 +y 2 



(73) 



Obviously, the introduction of proton quasiparticles af- 
fects only the proton part of the matrix element that de- 
scribes electromagnetic interaction. For the longitudinal 
part we have 

(2'|p p |2) = (u 2 'U2-V2'V 2 )5 <T2CT2( exp(z(p 2 /-p 2 )-r). (74) 

It is clear that the longitudinal component of the matrix 
element M| c differs from the same component M" orm in 
a nonsuperconducting (normal) case by 



A/| c = (u 2 'U 2 - v 2 'V 2 )M| 



(75) 



In the same way for the matrix element of the proton cur- 
rent (2'| J p \2), that describes the transverse interaction, 
we obtain 



(u 2 /u 2 -I- v 2 /v 2 ) M° 



(76) 



The differential transition probability for the decay 
channel has the form 



dec 
cp 



27r5(e 1 -e 2 -ev-e v ) ^ \ ,~2,2'\V\1) 

spins 



(77) 

Notice that the energy conserving delta function differs 
from that in the scattering process. However, momen- 
tum and spin selection restrictions are the same for scat- 
tering, decay and coalescence processes because of the 



adopted choice of initial and final particle momenta and 
spin quantum numbers. We obtain 



= (u 2 /v 2 + u 2 v 2 /) sign(cr 2 )M|] lorm , 
Ml cc = (u 2 ,v 2 -u 2 v 2 0sign(a 2 )Af7 rm . 

Similarly, for the coalescence channel 



(78) 



W. 



coal 



= 27T<5(£i + £ 2 - Ey 



■e 2 > 



J2 kl'1^1,2,-2') 



spins 



M 



coal 



coal 



(u 2 /v 2 + u 2 v 2 >) sign(cr 2 )M | ] lorm , 
(u 2 v 2 , - u 2 ,v 2 ) sign( ( r 2 )Mi orm . 



(79) 



(80) 



Our procedure is standard and tacitly ignores renormal- 
ization of proton vertices in superconducting matter. It 
is justified because we are interested is small energy and 
momentum transfers lo and q which are far from the 
transfers typical of collective superfluid excitations (see 

Refs. gi|[3fll and sec.nrsj. 

While integrating over quasiproton momenta in col- 
lision integrals we set d 3 p = m* p Tpp p dx dfl. It is in- 
structive to replace x 2 — > —x 2 for the decay case, and 
X2> — ► — X2> for coalescence. Such replacements change 
sign of quasiproton energy [e — fi p ) and coherence fac- 
tor v p because of our choice of these factors. Then all 
energy-dependent terms, except for coherence factors, be- 
come the same in all three collision integrals. Hence we 
can describe the interaction of electrons (or muons) and 
quasiprotons by a unified collision integral (|63[> . where 
W^(l,2|l',2') is substituted by Wg(l,2\l',2'), 



V^ P (l ) 2|l',2')=4(2 7 r) 4 %i + 
xS(p 1 +p 2 -p' 1 -P2)\M, 



e 2 - 

S |2 
cp I ' 



4) 



(81) 



\M, 



/fS\2 _ 



cp I 



+- 



M%%-x 2 )\ 2 



coal 



2 \ M cp 



(-4)| 2 . (82) 



Putting the squared matrix elements together, we have 



tp\\ = (1 - 4u 2 U 2 'V 2 V 2 /) i/3||, 

- (1 +4u 2 u 2 'V 2 v 2 /) <p±, 
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s 



(1 -vl-vj,)ip,n, 



(83) 



instead of tpn, tp±, and (p\n given by Eqs. (|TS|) — (TSTj) in 
the nonsuperconducting case. 

Therefore, the expressions for v ep or i> w in the pres- 
ence of proton superconductivity can be presented in the 
integral form similar to that in the nonsuperconducting 
case, but with two differences. First, coherence factors 
have to be introduced in the squared matrix element in 
accordance with Eq. ([83]) . Second, proper quasiproton 
energies (1731 . containing energy gaps, should be used in 
Fermi-Dirac distribution functions of quasiprotons. Thus 
in order to calculate v ep and v^p we should reconsider not 
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only the angular integration due to changes in polariza- 
tion functions but the energy integration as well. The 
collision frequencies in question contain the transverse 
and longitudinal components (|33[) . A careful examina- 
tion of the derivation of the energy integral in (|3Tj) shows 
that for the collisions via the exchange of transverse plas- 
mons one has 



cp 



(84) 



6C(3) 



dx 2 dx 2 
1 + cxp(z 2 ) 



(4 ~ Z 2)I4 - ^\FM}4 ~ z 2 \,y,r) 
[1 + exp(-4)][exp(4 - z 2 ) - 1] 

X(l + 4U2U2V2V2) 

(z' 2 + z 2 )\z' 2 + z 2 \F ± {\z' 2 + z 2 \,y,r) 



[1 + exp(z 2 )][exp(-z 2 

(1 - 4U2U 2 V2V 2 )| , 



z 2) 



where v^r corresponds to normal protons, Rp(y,r) de- 
scribes the effects of superconductivity, z 2 = (e 2 — 
H P )/k B T = yjx\ + y 2 and z' 2 = (e' 2 - fi p )/k B T = 

In the limit of strong superconductivity (y ^> 1) we 
have Rp — A±(r + 1) exp(— y), where 



A] 



2\3 



exp(ry?) - exp (77! 



0.3446. 



(85) 



This strong superconductivity exponentially reduces the 
ep and \xp collision rates. Nevertheless, the protons give 
the main contribution to the plasma polarization associ- 
ated with the exchange of transverse plasmons. In this 
way they remain vitally important for the plasma screen- 
ing in collisions involving leptons alone. 

The total transverse collision frequency of heat carriers 
c is 



(86) 



where v^r is the total collision frequency in nonsupercon- 
ducting matter and 

RL (V, r ) = — j— r [rRf (y, r) + R p L (y, r)] (87) 
r + 1 

describes the overall superconducting suppression of the 
total collision frequency via the exchange of transverse 
plasmons. In the limit of strong superconductivity, we 
havei? t L ot ( 2/ ,r)=4 7 r 2 r/[45C(3) 2 /]. 

We have calculated R^ t (y, r) on a dense grid of y and 



r and fitted the results by the formula 

2\ , 1 ~Pl 



i? t L ot (y,r) = pi exp (-0.1V 



P3V Z 



Pi 



0.4 



0.1 7r, 



P3 = 



(1-Pi) 



45C(3) 
47r 2 r 



which reproduces also the asymptote of R^ Qt (y,r) at 
y — > 00 and the obvious condition i?j^ t (0,r) = 1. The 
maximum fit error rj 5% takes place at y — 1.2 and r = 1. 

For the longitudinal components of the collision fre- 
quencies v\ v and vftp, the angular integration in (j25|) re- 
mains unchanged. Then 



4p=4 P RHy)> 

where v\ v is the nonsuperconducting value, while 

R\{y) = 



(89) 



00 00 

15 f f dx 2 dx 2 

4^7 J 1 





exp(^ 2 ) 

(4 - Z 2 ) 3 (l - 4u 2 U 2 V2V2) 

[1 + cxp(-z 2 )][exp(z 2 - z 2 ) - 1] 
(z 2 + z 2 ) 3 (l + 4u 2 u 2 v 2 v 2 ) 



[] +( . M)( -V)][ l . X p(-^- Z a)-l]J (90) 

A\\ exp(-y), 



is the superfluid reduction factor. 

In the limit of y 1 we obtain Rp{y) 
where 



15 
2^4 



d?7i / d?7 2 

(^-^) 3 (^+^2 2 ) 

exp(r?f ) - cxp^ 2 ) 



1.00. 



(91) 



Finally, in analogy with Eq. (|88p , we have numerically 
calculated R P {y) and fitted it by 

R\{y) = {l+(26.33y 2 + 0.376y 4 ) 
(-V(3.675) 2 + y 2 ) 

((1.673) 2 - V(l-673) 2 + y 2 ) - l] } 

(92) 



x exp 
+ 0.742 



exp 



exp ( (1.361) 2 - v/(1.361) 2 + y 2 



with the error < 1 %. 

The reduction factor iijj(y) was introduced by Gnedin 
and Yakovlev [17J [see their Eq. (43)]. In contrast to our 
Eq. ([90]) . their reduction factor had two drawbacks. First, 
it neglected coherence factors - the terms containing U2 
and V2 in (f9"0"|) . Second, the authors of Ref. [l7| attributed 

the same reduction Rp(y) to collisions via the exchange 
of transverse and longitudinal plasmons [whereas in fact 
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the reduction factors R^(y,r) and Rp(y) are different]. 
However, because superfiuid reduction of cp collisions for 
y 3> 1 is exponentially strong, these drawbacks do not 
change noticeably numerical values of the thermal con- 
ductivity (while neglecting the exchange of transverse 
plasmons makes results really inaccurate). 

Proton superconductivity significantly modifies the 
thermal conductivity of electrons and muons. In par- 
ticular, it violates the temperature independence of the 
conductivity (|40[) for nonsuperconducting matter. The 
asymptote (|87[1 indicates that in a strongly superconduct- 
ing matter (T -C T cp , |/> 1) the collision frequency 
becomes much smaller than in the nonsuperconducting 
case. Strong superconductivity enforces the temperature 
dependence v^r s oc T 2 that is formally the same as for 
the longitudinal collision frequency v\ in normal Fermi 
liquid. Nevertheless, this thermal conductivity remains 
smaller than the conductivity in the normal matter if 
it included artificially the exchange of longitudinal plas- 
mons alone. 

As in normal matter, the major contribution to the 
electron or muon thermal conductivity in the presence 
of proton superconductivity comes from the exchange of 
transverse plasmons. Using the asymptotes for the re- 
duction factors in the limit of strong superconductivity 
(y 3> 1) we obtain (in standard physical units) 



_5_ k B cp 2 Fc A 

24 ah 2 k B T p 2 F 



P 2 f p 



(93) 



Therefore, strong proton superconductivity formally re- 
stores the temperature dependence of the thermal con- 
ductivity Kefx oc 1/T typical for a normal Fermi liquid 
(by enforcing nondissipative plasma screening). Accord- 
ingly, we expect that in this case our simplest variational 
thermal conductivity is different from the exact conduc- 
tivity by the same factor C w 1.2 that is characteristic for 
Coulomb scattering via the exchange of longitudinal plas- 
mons [171 ]. These arguments support our expectations 
(Sec. Ill CI) that the overall uncertainties of our calcula- 
tions of k cm are ~ 20%. Notice that k~ m is proportional 
to the superfiuid gap A. 



IV. RESULTS AND DISCUSSION 
A. Equations of state 

Let us illustrate our results by calculating the thermal 
conductivity in a neutron star core composed of neutrons, 
protons, electrons and muons. We employ three repre- 
sentative equations of state (EOSs) which we denote as 
APR, PAL II, and PAL IV. The fractions of electrons 
x e = rie/rif, and muons = n^/rib versus density of 
matter pi4 — p/(10 g cm -3 ) for these EOSs are plot- 



ted in Fig. [H (n b 



being the number density 



0.01 



0.001 




of baryons). Let us recall that the lowest density in a 
neutron star core is p m ; n w 1.5 x 10 14 g cm -3 . 



Figure 2: (Color online) Fractions of electrons (thick lines) 
and muons (thin lines) versus density in a neutron star core 
for three selected EOSs — APR (dashed lines), PAL II (solid 
lines), and PAL IV (dash-dot lines). 



The APR EOS was constructed by Akmal, Pandhari- 
pande and Ravenhall [33j (their model Argonne V18 + 
Sv + UIX*). It is sufficiently stiff and gives the maxi- 
mum gravitational mass of stable neutron stars M max = 
2.2 Mq. The central density of the maximum-mass star 
is pmax = 2.8 x 10 15 g cm -3 . No hyperons appear for 
this EOS at p < p max . The threshold density for the 
appearance of muons is p M « 2.28 x 10 14 g cm -3 . 

The PAL II is a convenient phenomenological semi- 
analytical model EOS proposed by Prakash, Ainsworth 
and Lattimer [34|. The authors suggested several EOSs 
of such a type which differ by a value of the compres- 
sion modulus K of saturated symmetric nuclear matter 
and by the symmetry energy S of dense matter as a func- 
tion of baryon number density rib [described by a function 
F(u), where u = rib/rio and uq « 0.16 fm -3 is the baryon 
number density in saturated symmetric nuclear matter]. 
The compression modulus determines the stiffness of an 
EOS, while the symmetry energy regulates the fraction 
of protons (and hence of electrons and muons). The PAL 
II EOS corresponds to F(u) = 2u 2 /(u + 1) which gives 
a rather high fraction of protons at u > 1 and the muon 
threshold density p M = 2.647 x 10 14 g cm" 3 (Fig. [2J). 
Number densities of various particles, and hence K e/J , are 
insensitive to the value of Kq for PAL EOSs [34|]. Notice 
that by taking K o =120, 180, and 240 MeV (as suggested 
in Ref. [3J]), we would obtain three modifications of the 
PAL II EOS with different compressibility (from soft to 
stiff) which would give very different neutron star models 
(different mass-radius relations and maximum masses). 

The PAL IV EOS belongs to the same PAL family of 
EOSs [34] but corresponds to the function F(u) = u ' 7 
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Figure 3: (Color online) Electron-muon thermal conductivity 
K e(1 versus density in a nonsuperconducting neutron star core. 
Thick lines are for the APR (dashed line), PAL II (solid line) 
and PAL IV (dot-dashed line) EOSs at T = 10 8 K. Thin solid 
line is for the PAL II EOS at T = 10 9 K. 

suggested by Page and Applegate [35[ . It gives noticeably 
lower fractions of protons, muons and electrons at u > 1 
than those provided by the PAL II EOS (Fig. |2J. The 
muon creation threshold is p M = 2.51 x 10 14 g cm~ 3 . 
Taking again different K we get different modifications 
of the PAL IV EOS, from soft to stiff, giving different 
neutron star models. 

Therefore, the selected EOSs cover a large range of 
physical models of npep matter in neutron star cores. 



B. Thermal conductivity in nonsuperconducting 
cores 

Figure [3] shows the density dependence of the full 
electron-muon thermal conductivity K efi for the selected 
EOSs. The thick lines give at T = 10 8 K, a typi- 
cal internal temperature of a middle-aged neutron star 
(t ~ 10 4 — 10 5 yr) with no enhanced neutrino emission 
generated in its core (see, e.g., Ref. j5|). The conductiv- 
ity K efJ , increases with growing p because of the increasing 
amount of electrons and muons in dense matter [see Eq. 
(|40)) ] . The conductivity shows a noticeable kink at p = p^ 
associated with the appearance of muons. The conduc- 
tivities for the APR and PAL IV EOSs are close to each 
other because these two EOSs give very similar amounts 
of electrons and muons (Fig. [2]). At p > 4 x 10 14 g cm~ 3 , 
the conductivity for the PAL II EOS is noticeably higher 
owing to larger fractions of electrons and muons for the 
PAL II EOS at high densities. 

The thin solid line in Fig. [3] displays K efJl for the PAL 



Figure 4: (Color online) Electron (dashed line), muon (dot- 
dashed line) , electron-muon (solid line) thermal conductivities 
(calculated for T — 10 s K but almost temperature indepen- 
dent) and neutron conductivity (for T = 10 8 and 2 x 10 9 K, 
thick and thin dotted lines, respectively) versus density in a 
normal neutron star core for the PAL IV EOS. 



II EOS at a higher temperature T = 10 9 K typical of a 
young neutron star (a few months after its birth in the 
absence of enhanced neutrino emission; e.g., Ref. @). It 
is very close to the conductivity at T = 10 8 K confirm- 
ing our prediction (|4"0")) that n efl is almost independent 
of temperature in nonsuperconducting matter because of 
the importance of the Landau damping. The Landau 
damping owing to the exchange of transverse plasmons 
is the dominant mechanism of electron and muon scatter- 
ing; it fully regulates n efl in nonsuperconducting neutron 
star cores (Sec. IIIBj) . 

Figure 3] compares the thermal conductivity provided 
by different particles in a nonsuperfluid neutron star for 
the PAL IV EOS. We show the density dependence of 
the electron conductivity K e , the muon conductivity k m , 
and the electron-muon conductivity n efl calculated at 
T = 10 8 K with the notice that these conductivities 
are almost temperature independent. We also plot the 
neutron conductivity K n calculated using the results of 
Baiko, Haensel and Yakovlev for T = 10 s K (the 
thick dotted line) and 2 x 10 9 K (the thin dotted line). 
It is determined by collisions of neutrons with neutrons 
and protons via nuclear interactions and demonstrates a 
traditional Fermi liquid behavior (n n ot T -1 ). The ef- 
fective masses of neutrons and protons are set equal to 
0.7 of their bare masses, and the in-medium effects on 
the squared matrix elements of nucleon-nucleon scatter- 
ing are neglected. It is seen that the neutron contribution 
dominates throughout the neutron star core at T = 10 s K 
(contrary to the previous results [l3l. fl7|). However, we 
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Figure 5: (Color online) Temperature dependence of the 
electron-muon thermal conductivity K efl (thick lines) and the 
neutron conductivity k„ (thin lines) in a neutron star core 
with the APR EOS at p = 4 x 10 14 g cm -3 . Solid lines 
refer to nonsuperconducting matter, while dot-dashed and 
dashed lines are for matter with proton superconductivity 
(with T cp = 10 9 and 9 x 10 9 K, respectively). 



have n e p > K n in the core of a hot neutron star with 
T > 2 x 10 9 K. 



C. Thermal conductivity in superconducting 
matter 

Let us illustrate the effects of proton superconductivity 
on the electron-muon thermal conductivity K efl . Figure [5] 
shows K eM and k„ versus temperature for the APR EOS 
at p = 4 x 10 14 g cm" 3 assuming either T cp — (normal 
matter; solid lines), or T cp = 10 9 K (dot-dashed lines), or 
T cp = 9x 10 9 K (dashed lines). The neutron conductiv- 
ity is calculated using the results of Ref. [l5[ under the 
assumption that neutrons are normal. 

In normal matter K n dominates over K efl at T < 
2 x 10 9 K in agreement with the results of Sec. IIVBI We 
see that n efl is almost temperature independent because 
of the dominant contribution of the Landau damping [Eq. 
(|4T))) ]. Proton superconductivity reduces the electron and 
muon collision frequencies and increases K eLt as discussed 
in Sec. IIIIl it also slightly increases K n UM. Therefore, 
when the temperature drops below T cp and proton super- 
conductivity sets in, K e)Jt starts to grow up much quicker 
than K n and becomes comparable to or larger than K n . 
At T < T cp /3 this increase of K efl is well described by 
the asymptotic expression (jM)) . n efl oc T~ x . Proton su- 
perconductivity formally restores the Fermi liquid be- 
havior of (Sec. IIIIl) . However, this K efl is smaller 



than the electron-muon conductivity which is calculated 
in Refs. [H, Nevertheless, because K efl in supercon- 
ducting matter is directly proportional to the proton gap 
A oc T cp , the enhancement of K eM by proton supercon- 
ductivity makes n efi quite large. According to Fig. [5l for 
T cp — 10 9 K and T < 3 x 10 8 K we have K efl much larger 
than in normal matter although a factor of ~ 3 lower than 
K n . For T cp — 3 x 10 9 K (not shown in the figure) we 
would have re eM ~ K n at any T displayed. For a stronger 
superconductivity with T cp = 9 x 10 9 K, as seen from Fig. 
K e fi dominates over re„ at any T. Therefore, we obtain 
K efl > K n for T > 2 x 10 9 K in normal matter and for 
any T in superconducting matter with T cp >3x 10 9 K. 

V. CONCLUSIONS 

We have calculated the thermal conductivity of elec- 
trons and muons K e/J , in the cores of neutron stars com- 
posed of neutrons, protons, electrons, and muons taking 
into account possible superconductivity of protons. This 
thermal conductivity is determined by electromagnetic 
interactions of electrons and muons with all charged par- 
ticles. Our results are presented in the form of simple 
analytic expressions and fitting formulae valid for any 
equation of state of dense matter in a neutron star core 
and for any ratio T/T cp between the temperature T and 
proton critical temperature T cp in superconducting mat- 
ter (T < T cp ). A generalization of our results to hyper- 
onic matter (including the case of superfluid hyperons) 
is straightforward (analogous to that in Ref. p7|). 

Our main conclusions are as follows. 

1. The main contribution to K e/i comes from electro- 
magnetic interactions of charged particles via the 
exchange of transverse plasmons. This contribu- 
tion has been neglected in all previous calculations 
of K e p in neutron star cores (although it has been 
included in the calculation of the thermal conduc- 
tivity of quarks in quark matter [l8| and the ther- 
mal conductivity of electrons in neutron star crusts 

My 

2. For normal (nonsuperconducting) protons, n ell is 
determined by electromagnetic interactions of elec- 
trons and muons with all charged particles (e, p, 
p) via the Landau damping of transverse plasmons 
[Eq. PD)) ]. This thermal conductivity is tempera- 
ture independent (contrary to the traditional Fermi 
liquid behavior k oc T~ x ). 

3. The conductivity n e p is mainly determined by the 
symmetry energy of dense matter (rather than by 
the stiffness of the equation of state) and increases 
with the growth of the symmetry energy. 

4. At temperatures T ~ 10 8 K, the conductivity K efl 
in a normal neutron star core is smaller than the 
conductivity of neutrons K n , but at T > 2 x 10 9 K 
we obtain K efi > n n . 
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5. After the onset of proton superconductivity (at 
T < T cp ) remains to be determined by electro- 
magnetic interactions of charged particles via the 
exchange of transverse plasmons but the character 
of these interactions is different. First, strong pro- 
ton superconductivity greatly suppresses the colli- 
sions of electrons and muons with protons. Sec- 
ond, this superconductivity modifies the exchange 
of transverse plasmons (from the dissipative Lan- 
dau damping to the nondissipative regime). As 
a result, strong superconductivity (T < T cp /3) 
increases K efJ _, restores its Fermi liquid behavior, 
and enforces n efl to scale as K eM cx T cp /T. For 
T cp >3x 10 9 K, we obtain K e ^ > n n at any T. 

Our results are based on previous calculations of 
Gnedin and Yakovlev [TtJ • We have extended them by 
properly including the charged particle interactions via 
the exchange of transverse plasmons and by performing 
exact calculation of the rate of electron/muon scattering 
by proton quasiparticles in the presence of proton super- 
conductivity. 

We expect that our results give reliable values of K efl 
in neutron star cores. In the absence of proton supercon- 
ductivity, calculation of K eM constitutes a well defined 
problem, which we solve explicitly (under the only one 
and well justified assumption of weak plasma screening). 
Because the solution involves only electromagnetic inter- 
actions, it is universal. It does not depend on a strong 
interaction model of dense matter and, hence, on a spe- 
cific equation of state in a neutron star core. In the pres- 
ence of proton superconductivity, we use the polarization 
functions of a proton gas derived in the BCS framework. 
We argue (Sec. IIII A|) that it can be a good approxima- 
tion. Even if future studies of realistic nucleon Fermi 
liquid give more elaborated polarization functions, our 
solution will serve as a useful basic approximation. 

In any case our results do not solve the thermal conduc- 
tion problem in neutron star cores. The main problem is 
posed by the conductivity of neutrons and protons (and 
other baryons if available) mediated by strong (nuclear) 
interactions, especially in the presence of baryon super- 
fluidity. The conductivities of neutrons and protons, n n 
and n pi in normal matter have been estimated by Flow- 
ers and Itoh [H, [3 (with the natural result that ftp is 
much smaller than K n because of smaller amount of pro- 
tons in neutron star cores). More accurate estimations 



of K n were done later by Baiko et al. [15[ who estimated 
also the diffusive thermal conductivity k„ in the pres- 
ence of neutron and proton superfluidity. However, the 
appearance of neutron superfluidity can trigger a specific 
and efficient nondiffusive heat conduction via convective 
counterflow. This effect is well known from laboratory 
experiments with superfluid 4 He (e.g., Ref. [36|): it is so 
efficient that immediately dissolves any temperature gra- 
dients in superfluid 4 He. Analogous effects in superfluid 
neutron star cores have been mentioned in the astrophys- 
ical literature (e.g., Ref. [l3j]) but have not been explored. 
Their careful examination would be desirable. 

It would also be important to study the effects of strong 
magnetic fields which can greatly modify thermal con- 
ductivity in neutron star cores. Strong Larmor rotation 
of charged particles (electrons, muons, protons) about 
magnetic field lines can greatly suppress thermal con- 
duction across the magnetic field. Similar effects in elec- 
tric conductivity have been studied in a number of works 
(see, e.g., [37j and references therein) but they are almost 
unexplored for the thermal conductivity. The magnetic 
field modifies plasmon modes in dense matter and affects 
plasma polarization properties, particularly, the Landau 
damping. The magnetic field effects can be especially 
complicated in the presence of baryon superfluidity. 

All in all, further serious efforts are required to solve 
the heat conduction problem in neutron star cores. First 
of all, its solution is needed to model cooling of neutron 
stars [E [f| (especially in the first 100 years of their life 
[ToL [llf), thermal states of neutron stars in soft X-ray 
transients [1, 0, H, Q and thermal relaxation of pulsars 
after glitches [12| . We hope that our results give a reliable 
contribution to the thermal conduction problem, and we 
expect to study other effects of the Landau damping on 
kinetic properties of neutron star cores (particularly, on 
the shear viscosity) in subsequent publications. 
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